Ab initio potential curves for the X 2 £+, A 2 IT M and B 2 Z+ states of Ca^ 

Sandipan Banerjee*, John A. Montgomery, Jr., Jason N. Byrd, H. Harvey Michels, Robin Cote 

Department of Physics - University of Connecticut, Starrs, CT 06269-3046, USA. 



(N 



Abstract 



We report ab initio calculations of the X 2 £+, A 2 Tl u and B 2 E+ states of the Ca} dimer. All electron CAS+MRCI calculations 
are performed for the X 2 £+ and B 2 £+ states, while valence CAS+MRCI calculations using an effective core potential are used 
^ to describe the A 2 II„ state. A double well is found in the B 2 E+ state. Spectroscopic constants, vibrational levels, transition 
^ moments and radiative lifetimes are calculated for the most abundant isotope of calcium ( 40 Ca). The static dipole and quadrupole 
H~5 polarizabilities, and the leading order van der Waals coefficients are also calculated for all three states. 
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Introduction 

The presence of near degeneracies in the constituent atoms of 
diatomic molecules can lead to a rich structure in the resulting 
interaction potentials. In our recent work on the Be} dimer (Q]], 
we showed that the nearly degenerate 2s-2p state of beryllium 
leads to a complex set of low lying molecular curves, including 
a double minima in the lowest 2 S^ state. In this Letter, we 
present new computational results on the Ca} dimer, and 
demonstrate that this system also exibits a rich manifold of 
low-lying molecular states. Both Be and Ca represent group 
2A elements in the periodic table, with their valence electronic 
structures (2s) 2 and (4s) 2 , respectively. Thus, some similarities 
between Be} and Ca} may be expected. The calculations 
presented here should help guide experimental efforts on cold 
molecular ions. 

The last few years have seen significant interest in ultracold 
atom-ion scattering Si in the atomic, molecular and optical 
physics community. The experimental realization of Bose- 
Einstein condensation (BEC) has led to numerous applications 
involving charged atomic and molecular species. The cooling 
and trapping [4] of charged gases at sub-kelvin (ultracold) 
temperatures is a topic of growing interest. The phenomena of 
charge transport like resonant charge transfer Q5J and charge 
mobility J(J at ultracold temperatures have also been studied 
in detail. Other emerging fields of interest include ultracold 
plasmas 10], ultracold Rydberg gases and systems involving 
ions in a BEC fflfioll. 

In this work, we begin by describing the methods used in our 
calculations, and follow with a discussion of the results, which 



include the potential curves of the X 2 X,}, B 2 £+ and A 2 H U 
states and their spectroscopic constants. We also calculate elec- 
tric dipole transition moments for the X 2 E} <-> B 2 E+ and the 
B 2 l,g <-» A 2 n„ transitions. Bound vibrational levels are com- 
puted for all the states along with Franck-Condon overlaps and 
radiative lifetimes for the most abundant calcium isotope ( 40 Ca 
96.94%). We conclude with an analysis of long range behavior, 
calculation of static atomic dipole and quadrupole polarizabil- 
ities and determination of the van der Waals dispersion coeffi- 
cient Ce- 



Methods 

We can express the total energy of the Ca} dimer at any inter- 
atomic separation R as 
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states in Ca}, 



For calculation of the ground X 2 X+ and B 2 SJ 
the valence contribution to the total energy is calculated by 
a multi-reference configuration interaction (MRCI) method 
using a 18 orbital complete active space (CAS) wavefunction as 
a reference. The active space was chosen to include molecular 
counterparts of nearly degenerate 4s, 4p and 3d orbitals of Ca. 
The state-averaged CAS includes all doublet states correlated 
to Ca + ( 2 D) + Ca( r S) and Ca + ( 2 S) + Ca('S) asymptotes with 
equal weights. We have used the augmented correlation 
consistent polarized valence quintuple zeta (aug-cc-pV5Z) 
basis set of Peterson 111 ll Il2ll . In order to assess the quality 
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of MRCI, we do a comparison with a full CI calculation with 
aug-cc-pVTZ basis, and find out that the difference in total 
energy at the equilibrium separation of 7.3 bohrs for the X 2 Z+ 
state of CaJ is 4.5 microhartrees. At large separation (1000 
bohrs), this difference further reduces to 1.5 microhartrees. 
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Figure 1: [COLOR ONLINE] Ab initio X 2 Z,+ (in black), B 2 X+ (in red) and 
A 2 n„ (in blue) states of Ca+. Dashed lines show calculation with a pseudopo- 
tential, while solid lines show the results of an all electron correlated calcula- 
tion. For the A 2 n„ state, results of all electron MRCI calculations are shown 
in green •. Note that ao is the Bohr radius (atomic unit of length). 



The second term in Eq. (Q}, the correction from the core- 
valence contribution is estimated by, 
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The core-valence correction AE core - va \ ence is the difference 
of energies from a valence only (Ar core) and a restricted 
inner valence (RIV, Ne core) CCSDT calculation. For this 
purpose we have used the correlation consistent polarized 
weighted core-valence triple zeta (cc-pwCVTZ) basis set of 
Koput and Peterson fiUl . To assess the convergence with basis 
set of the calculated core-valence contribution, we performed 
single point calculations at the equilibrium bond separation 
for the ground state with a larger basis set (cc-pwCVQZ). The 
effect of increasing the basis set from TZ to QZ changed the 
core-valence energy by ~ 8 crrT 1 at the R e of the B 2 E+ state. 
These results indicate that the core-valence contribution to the 
total energy is adequately converged with the TZ basis sets. 

The last correction term A£ sca / ar _ re toj V , J( ,- c is the contribution 
from relativistic effects, which for a heavy atom like Ca is sig- 
nificant. This can be expressed as, 
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We have used the Douglas-Kroll version of the cc-pwCVTZ 
basis set from Kirk Peterson (cc-pwCVTZ-DK), and performed 
CCSDT calculations to estimate this correction. The magnitude 
of scalar relativistic correction at the equilibrium bond distance 
for the B 2 E+ state of the CaJ is ~ 180 crrT 1 . For Ca+, the 



valence electron space contains only 3 electrons, thus the 
valence CCSDT is equivalent to full CI. 

The calculation of the A 2 LI„ state, correlating to the Ca + 3d 
atomic level, is complicated by the near degeneracy with the 
Ca 4s4p atomic level. The second excited 2 II„ state comes 
from an atomic asymptote of Ca 4s4p and Ca + 4s, which lies 
~ 1500 cirT 1 above the A 2 lL. asymptote. We find, however, 
that valence CAS+MRCI calculations incorrectly predict the 
Ca 4s4p and Ca + 4s asymptote to lie below the Ca 4s 2 and 
Ca + 3d asymptote. The correct ordering of the atomic energy 
levels is obtained when core-valence correlation including 
double excitations of the inner valence electrons are included 
in the correlation treatment using the cc-pwCVQZ (or better) 
basis set. We expect that a balanced description of valence and 
core-valence interactions in the A 2 n„ state would be obtained 
from a CAS(19,26)+MRCI calculation that includes molecular 
orbitals arising from the atomic 3s, 3p, 3d, 4s and 4p orbitals. 
This is a much more demanding calculation than those required 
for the X 2 S^ and B 2 2J states that correlate to ground states 
atoms. It was found that a smaller CAS(19,21)+MRCI+Q/cc- 
pwCVQZ (MRCI plus Davidson correction) calculation 
correlating the 3s3/;> inner valence electrons with a 4s3d + 4p x 
valence reference was sufficient to obtain the correct ordering 
of the first and second 2 Yl u molecular states. Extending this 
calculation to include the entire 4p reference space using 
the cc-pwCV5Z was attempted but was too computationally 
demanding for our available resources. An alternative ap- 
proach to the multi-reference all electron calculation is to 
replace the argon core of the Ca atoms with an effecive-core 
potential (ECP) where the effects of core-valence correlation 
are including using a core polarization potential (CPP) I113I1 . 
This method was used with great success by Czuchaj et al 
for Ca2 ground and excited states 111 411 . We have performed 
comparisons between the all electron calculations for the first 
two 2 S states as described above and the valence 4s4p3d 
space MRCI using the ECP+CPP and basis set of Czuchaj et 
al 01411 , The agreement was found to be satisfactory for the 
case of the lowest S states as seen in Fig. [TJ Additionally 
we have compared the A 2 LI„ state calculated using the same 
ECP+CPP method to the core-valence MRCI+Q/cc-pwCVQZ 
calculation using the 3s3p4s3d + 4p x space discussed above. 
These two calculations agree very well, as demonstrated by 
Fig. [TJ Because of the good agreement with the all electron 
calculations and the computational limitations in performing 
multi-reference core-valence correlation calculation, we have 
used the ECP+CPP method to calculate the A 2 n„ state in this 
paper. We note in passing that the use of the CPP is essential; 
without it one does not obtain the correct ordering of the 
excited atomic asymptotes. 

All the potential curves are also corrected for basis-set super- 
position error (BSSE) using the standard counterpoise tech- 
nique of Boys and Bernardi [15]. The BSSE was negligible 
(~ 2 - 4 cirT 1 ) at the potential minima for the different curves. 
The MRCI valence calculations were done using the M0LPR0 
2010.1 electronic structure program lll6l . The core- valence 
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Figure 2: [COLOR ONLINE] Calculated ab initio potential curves of CaJ. 
The inset is a magnification of the double-well nature in the B 2 Z+ state (in 
red); lowest vibrational levels in the inner well are shown in blue thick lines 
and outer well in red thin lines, for 40 Ca. 



Table 1 : Calculated spectroscopic constants of Ca+ 



State 


r,(A) 


B e (cm -1 ) 


w e (cm ') 


cu e x e (cm ') 


D e (cm -1 ) 




3.844 


0.056 


127.829 


0.071 


9440 


Previous f 141 


3.773 




132.300 




9817 


Previous [231 


3.995 


0.053 


119.000 




8388 


B 2 I+ (Inner) 


4.719 


0.037 


41.593 


0.561 


284 


B 2 Z+ (Outer) 


8.665 


0.011 


8.549 


3.839 


137 


A 2 n„ 


3.303 


0.077 


194.195 


0.370 


14746 



CCSDT calculations were carried out using CF0UR (coupled- 
cluster techniques for computational chemistry) program lfl7ll . 
The scalar relativistic corrections were done at the CCSDT level 
of theory using the MRCC (multi-reference coupled cluster) pro- 
gram [18] of M. Kallay. All of the programs were running on a 
Linux workstation. All calculations employed restricted open- 
shell (ROHF) reference wavefunctions. Le Roy's LEVEL pro- 
gram U9A has been used to calculate the bound vibrational lev- 
els, Franck-Condon factors and radiative lifetimes, discussed in 
the following section. 



Results and Discussions 

Potential Curves and Spectroscopic Constants 

Fig. |2] shows the ab initio potential curves for the X 2 E+ and 
B 2 X+ states of Ca^ dimer. The calculated potential energy 
curves are corrected for the effects of basis set superposition 
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Figure 3: [COLOR ONLINE] (a) Calculated radiative lifetimes of bound levels 
of 40 CaJ in the B 2 X+ state, on a log-scale. The shorter lifetimes (blue •) 
correspond to bound levels localized in the inner well, the longer lifetimes (red 
o) to levels localized in outer well, and the increasingly longer lifetimes (green 
♦) to levels spread over both wells, (b) show examples of each cases: v=0 and 
7 in the inner well, v=6 in the outer well, and v=40 in both, (c) depicts |i//-| 2 of 
v=6 on a log-scale; showing that the amplitude in the inner well is extremely 
small. The reverse is true for v=7 and so on. 



error by the counterpoise method of Boys and Bernardi[15!|. 
FigQ] shows the A 2 Yl u state. We have used a standard Dunham 
analysis i20ll to calculate the spectroscopic constants (Table[TJ. 
We calculate bound vibrational levels for the X 2 £+ B 2 ZJ and 
A 2 IT„ state for the ^CaJ dimer. 

Unfortunately there are no experimental spectroscopic data 
for the ground or excited states of the Ca^ dimer. There are, 
however, so me p revious theoretical studies of Ca2 lfl4L Elll 
and Ca} l22l |23J, and a comparison to the results for the 
X 2 £+ state is listed in Table \T\ No spectroscopic constants 
have been reported for the most recent calculation by Sullivan 
et al. 022I1 . Another approach for calculating interaction 
energies in alkaline earth elements, is using symmetry-adapted 
perturbation theory (SAPT), which has been demonstrated 
earlier by Patkowski et aZ II24ll . The B 2 S^ state has a double 
well similar to that found in our recent Bet calculations |l|]. 
Both of these wells support bound vibrational states. This 
double-well nature of the B 2 2^ state is most likely caused by 
perturbations from an excited 2 E^ state. 

We have calculated radiative lifetimes (see PigOJ for bound vi- 
brational levels in both the inner and outer wells of the B 2 2^ 
state. Since these wells are separated by a large barrier, the 
wavefunction of the lower vibrational levels can be strongly lo- 
calized in either wells. The localization of the vibrational wave- 
functions can be attributed to the asymmetry of the double well 
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Figure 4: [COLOR ONLINE] Energies of bound levels of B 2 Z+ state in 40 CaJ 
using the same convention as in Figf5] The inset magnifies the difference in 
slopes of levels localized in the inner well (shown in blue) from the ones in the 
outer well (shown in red). 

(see FigO and disappears for levels above the barrier. The be- 
havior of the radiative lifetimes exhibits 3 distinct regimes, one 
where the wavefunction is mainly localized in the inner well 
(in blue), one where it is localized in the outer well (in red), 
and the last being the region (in green) in which the wavefunc- 
tion spreads over both wells, resulting in poor Franck-Condon 
overlap with the ground X 2 £+ state and hence longer lifetimes. 
The inset shows the square of the amplitude of wavefunction in 
v=6, of B 2 £J state in a logarithmic plot as a demonstration that 
the amplitude is negligible inside the inner well but still finite, 
preserving the correct number of nodes for that level. Fig@] 
shows a plot of the energies of all bound vibrational levels in 
the B 2 Z^ state of 40 Ca2 . The localization effect of wavefunc- 
tions discussed above is also exhibited in this plot; the density 
of levels in the more extended outer well is larger than in the in- 
ner well, leading to different energy slopes. The inset of Figj4] 
exemplifies this point. 

Electronic dipole transition moments 

For homonuclear molecules like CaJ, there is no permanent 
dipole moment. However there are transitions between differ- 
ent electronic states which are dipole allowed. We calculate 
two transition moments, one of them couples the ground 
X 2 Z+ and B 2 2^ states, and the other one couples the B 2 £+ 
state to the excited A 2 n„ state. To compute these transition 
moments, we use a complete active space self consistent 
field (CASSCF) wavefunction as a reference for performing 
multi-reference configuration interaction (MRCI) calculations. 
The core-valence contribution to the electronic transition 
moment is found to be negligible and hence omitted in the 
present calculations. The calculation of the transition moment 
coupling the B 2 £+ state to the excited A 2 n„ state was done 
with the ECP+CPP valence CAS+MRCI method. 

The electronic dipole transition moment (in atomic units) is 
given by, 




(i 12 (R) = (2\z\l), 



(4) 



where |1) and |2) are the electronic wave functions correspond- 
ing to the pair of states X 2 £+ <-» B 2 S^ or B 2 £+ <-> A 2 n„, small. Also we have truncated the series at powers of R b , not 



Figure 5: [COLOR ONLINE] Computed electronic dipole transition moment, 
fiBA coupling the B 2 Z+ to the A 2 n„ state shown in red (a) and fixB coupling 
the X 2 ZJ to the B 2 Z+ state shown in black (b). The dashed line R/2 in (b), 
corresponds to the classical dipole behavior. 



when the two Ca nuclei are separated by the distance R. 

Fig. |5ja) shows the electronic transition dipole moment 
coupling the B and the A II M state. The transition moment 
goes to zero asymptotically. Fig. Qb) shows a plot of the 
computed electronic dipole transition moment between the 
B 2 2g and the X 2 E+ ground states of Ca^. The transition mo- 
ment hxb coupling the X 2 £+ and B 2 £+ states asymptotically 
follows the classical dipole behavior, fixB ~ R/2 11251 12611 : we 
observe this behavior in the calculated curve of Fig. |5]b). We 
note that although the transition moment grows linearly with 
R, the probability of spontaneous transition will tend to zero 
since it is proportional to v XB , which vanishes exponentially as 
R — > oo. 



Polarizabilities and dispersion coefficients 

For large internuclear separations, the long-range form of the 
intermolecular potential can be written as 

V«CR) = V--2S. ( 5 ) 



which, for the molecular ion CaJ , can be approximated by 

V«W~V--§-§. (6) 

where Voo is the energy of the atomic asymptote; Cn—ad/2, ad 
is the static dipole polarizability, C(, = (a q /2 + c&), a q is the 
quadrupole polarizability and C(, the Van der Waals dispersion 
coefficient. In the expression for long range energy we have 
ignored the exchange energy contribution E exc f, which is very 

-6 
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Table 2: The static atomic dipole, quadrapole polarizabilities and dispersion coefficient for X 2 X+, B 2 X+ and A 2 n„ states of CaJ. All values are in atomic units. 
The square brackets indicate powers of ten. 



Molecular State 


Dipole 
polarizability («</) 


Quadrupole 
polarizability {a q ) 


c 4 

(=a d /2) 


Dispersion 
coefficient (c(,) 


c 6 

(=ff 9 /2+c 6 ) 


X a ZJ,B a 2£ 
Previous [271 
Previous [281 


1.606 [2] 
1.571 [2] 


3.073[3] 
3.081[3] 


8.032[1] 


1.081 [3] 
1.085[3] 


2.618[3] 


A 2 n„ 

Previous [301 


1.606[2] 


3.073[3] 


8.032[1] 


4.950[2] 


2.031[3] 



including contributions from R 8 and R 10 order coefficients. 

We have performed finite-field CCSD(T) calculations with the 
aug-cc-pV5Z basis set using M0LPR0 to estimate the values of 
the static atomic dipole and quadrupole polarizabilities. We ob- 
tain a c { = 160.64 a.u. and a q = 3073.39 a.u which are both 
in good agreement with a previous result jzzl of 157.1 a.u. 
and 3081 a.u. respectively. The values of static dipole and 
quadrupole polarizability do not change for the excited A 2 IT„ 
state, since it comes from an atomic asymptote in which the Ca 
ion is in an excited 2 D state whereas the Ca atom is in ground l S 
state (see FigJTJ. Using a numerical fit, as described by Baner- 
jee et al. HI, we were able to extract the value of the dispersion 
coefficient C6 for all the states. The value of dispersion coeffi- 
cient for ground states are in good agreement with unpublished 
results of Mitroy II28||. which are done by the methods used by 
Mi troy and ZhangUl 0] for Ca and Ca + . Table lists the 
values of polarizabilities and dispersion coefficient for X 2 £+, 
B 2 



i E+ and A 2 n„ states of the Ca^ dimer. 



Concluding Remarks 

Ab initio calculations have been performed on the X 2 £+, B 2 X+ 
and A 2 I1„ states of the CaJ dimer. The calculations were 
computationally challenging as well as expensive because of 
the near degeneracy of the 4s-4p-3d orbitals in Ca. Since the 
B 2 2g state has a double well, one of them being near 16.5 
ao, it was necessary to include diffuse functions in the basis 
sets to describe the well accurately. For calculating the ground 
X 2 T.+ and B 2 2^ states, large augmented correlation consistent 
basis sets of Koput and Peterson Hill were thus chosen and 
the results were also corrected for basis-set superposition 
error. We have also corrected our valence only MRCI results 
for core-valence and scalar relativistic effects using CCSDT 
calculations with both restricted inner valence and frozen core 
using Peterson's cc-pwCVTZ basis set. The A 2 JT H state was 
calculated with an ECP + CPP approach similar to the previous 
calculations of Czuchaj et al. 1 1411 . Numerical values of the 
calculated potential energies are available from the authors 
upon request. 

Due to lack of experimental data, we were unable to compare 
our theoretical values for dissociation energies or spectroscopic 



constants. However there are some previous theoretical results 
for spectroscopic constants for the X 2 Z+ ground state of Ca£ 
lfl4l l23ll . These values compare well with our calculated results 
(see Table[T). 

The calculation of radiative lifetimes of 40 Ca2 in the B 2 £+ 
state, reflects the wavefunction localization in either of the dou- 
ble wells. This causes the levels in the inner well to have a 
shorter lifetime (~ yus), whereas the ones in the outer well have 
much longer lifetimes (~ ms). These vibrational levels in the 
B 2 Z^ state should generate interest for experiments in ultracold 
atomic and molecular physics, where these long lived molecu- 
lar ions could be observed. We believe there are new prospects 
in both theory and experiments for atom-ion collisions, res- 
onant charge transfer and quantum information storage using 
Caj molecular ions. 
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